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We present evidence for the cooling of normal metal phonons, in addition to the well-known 
electron cooling, by electron tunneling in a Superconductor - Normal metal - Superconductor tunnel 
junction. The normal metal electron temperature is extracted by comparing the device current- 
voltage characteristics to the theoretical prediction. We use a quantitative model for the heat 
transfer that includes the electron-phonon coupling in the normal metal and the Kapitza resistance 
between the substrate and the metal. It gives a very good fit to the data and enables us to extract 
an effective phonon temperature in the normal metal. 



In a Normal metal - Insulator - Superconductor (N-I- 
S) tunnel junction, the presence of the superconducting 
energy gap induces the selective tunneling of high energy 
quasi-particles out of the normal metal [2, 0] • This phe- 
nomenon generates a heat current from the normal metal 
to the superconductor that is maximum at a voltage bias 
just below the superconductor energy gap. Only sin- 
gle electron tunneling contributes to the thermal current, 
while Andreev reflection does not. A double tunnel junc- 
tion (S-I-N-I-S) geometry optimizes the cooling power so 
that the electronic temperature of the normal metal can 
be lowered from 300 mK down to about 100 mK [3]. This 
significant temperature reduction raises hope to use N-I- 
S junctions for on-chip cooling [4] of nano-sized systems 
like high-sensitivity detectors and quantum devices. 

A complete investigation of thermal effects in such hy- 
brid nanostructures at very low temperature requires a 
precise determination of both the electron and phonon 
temperatures. A current-biased N-I-S junction on the 
normal metal can be used as an electron thermometer 
by measuring the voltage at a constant bias current [5[. 
However, the phonon temperature is very difficult to 
measure directly and is usually assumed to be that of 
the substrate. As we will detail below, this assumption 
hinders a full characterization of the thermal effects oc- 
curring in the device. Specifically, compared to previous 
studies @,0], the electrons are cooled here, which in turn 
can cool the phonons. The importance of this effect de- 
pends on the various energy relaxation processes, i.e. on 
the various thermal couplings between electrons, normal 
metal phonons and substrate phonons at the nanometer 
scale. 

In this Letter, we report on a precise investigation of 
the current- voltage characteristics of S-I-N-I-S junctions 
where the coupling between the normal metal phonons 
and the substrate is weak. The electronic temperature 
is obtained directly from the current-voltage character- 
istic. We obtain a quantitative description of the data 
with a heat transfer model that includes the electron- 




FIG. 1: Scanning electron micrograph of a cooler similar to 
Sample 1 featuring two large area Al-AlO^-Cu junctions in 
series. The inset shows the complete device with 3 probes on 
the superconducting side of the cooler junction. 



phonon scattering in the normal metal and the Kapitza 
resistance between the metal and the substrate. This 
analysis provides an effective phonon temperature in the 
normal metal which until now has never been accessed in 
such devices. 

In the following, we will discuss data from two typ- 
ical samples with a slightly different geometry. Fig. 1 
shows one of our S-I-N-I-S coolers which were fabricated 
using electron beam lithography, two- angle shadow evap- 
oration and lift-off on oxidized Si wafers. The two 40 
nm thick and 1.5 fim wide superconducting Al electrodes 
were in-situ oxidized in 0.2 mbar of O2 for 3 min before 
the deposition of the d = 50 nm thick, L =5 (Sample 
1) or 4 (Sample 2) fim long and 0.3 fim wide normal 
metal Cu electrode. The elastic mean free path in Cu 
is 33 nm. The two symmetric junctions of dimensions 
1.5 x 0.3 jim 2 give a total normal-state resistance in the 
range 2-3 kfi. The Cu electrode has been made rather 
short and thick. The thermal coupling between electrons 
and phonons is therefore increased relative to the cou- 
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FIG. 2: (Color online) Differential conductance data from 
Sample 1 measured at a cryostat temperature of 278 mK. 
Top curves: data (full blue line) from the middle external 
probe fitted with the isotherm (dotted red line) at 300 mK 
with R n = 5.7 kQ, A = 0.22 meV. Bottom curves: cooler 
junction data (full black line) compared with the isotherm 
model (dotted red line) at 300 mK and the electron-phonon 
model (dotted green line) with R n = 2.8 kQ, A = 0.21 meV, 
KA — 75 pW.K~ A . The voltage is normalized to A (top 
curves) and 2 A (bottom curves). 



pling to the substrate. The overlapping extra copper 
shadows outside the N-I-S junctions are assumed to act 
as quasiparticle traps for a better thermalization of the 
superconducting reservoirs [8]. In addition to the two 
cooling junctions, we added three Cu tunnel probes of 
area 0.3x0.3 fim 2 on one Al electrode (see Fig. [T] inset). 

We measured the current-voltage (I-V) characteristics 
of several S-I- N-I-S coolers in a He 3 cryostat. The differ- 
ential conductance dl/dV of every junction was obtained 
numerically. For sample 1, Fig. [2] shows, on a logarithmic 
scale, the differential conductance at a 278 mK cryostat 
temperature of both the middle external probe junction 
and the cooler double junction. The tunnel current of a 
N-I-S junction can be written as: 

1 f°° 

I(V) = — / n a (E)\f N (E-eV/2)-f N (E+eV/2)]dE, 

(i) 

where R n is its normal state resistance, /n is the elec- 
tron energy distribution function in the normal metal and 
n s (E) =| E | /VE 2 — A 2 is the normalized BCS density 
of states in the superconductor. Thus the tunnel current 
in a N-I-S junction does not depend on the supercon- 
ductor quasiparticle distribution. The probe data is well 
fitted by Eq. (1) in the hypothesis of a distribution func- 
tion /n at thermal equilibrium. The fit temperature of 
300 mK is only slightly higher than the measured cryo- 



stat temperature. The conductance at very small bias 
diverges slightly from the model while its linear slope 
(in log scale) versus the voltage in the sub-gap region 
agrees very well with the prediction. This good agree- 
ment demonstrates that the sample tunnel junctions are 
of good quality. 

Whereas the normal metal electrode of the probe junc- 
tion is strongly thermalized to the substrate temperature 
because of its large area, the central island of the cooler 
is small and thus weakly coupled to it. As can be seen 
in Fig. [2j the cooler junction indeed shows a strikingly 
different behavior that could not be fitted by Eq. (1) as- 
suming thermal equilibrium at any given temperature. 
The experimental curve and the calculated isotherm at 
300 mK coincide at zero bias, which confirms that the 
electron temperature is approximately equal to the cryo- 
stat temperature in the absence of current. In the subgap 
region, the differential conductance of the cooler is always 
below the isotherm. This behavior exemplifies clearly the 
cooling of the electronic population of the normal metal 
part of the device. 

The phase-coherence time in Cu was measured by 
a weak-localization experiment in a wire fabricated with 
the same material. This time gives an estimate for the 
electron-electron scattering time. We found = 0.19 ns 
at 275 mK, which is much smaller than the mean escape 
time of an electron out of the metal island estimated to be 
about 100 ns. We conclude that the electrons in the nor- 
mal metal island are at thermal quasi-equilibrium and 
that an electronic temperature T e can be well defined. 
We note that this regime is more easily met at the rel- 
atively high temperatures on which we concentrate here 
than at lower temperatures where deviations from quasi- 
equilibrium are expected to be prominent [9]. 

In the following, we assume that the two refrigerating 
N-I-S junctions are symmetric so that there is an equal 
voltage drop across each junction. Moreover, we assume 
the absence of leakage channels through the junction and 
of sub-gap energy states in the superconductor [10]. The 
Andreev reflection contribution to the subgap current 
[Tlj is here negligeable. We extract the temperature 
T e (y) in the sub-gap region by superimposing the ex- 
perimental I-V curve on a series of isotherm curves ob- 
tained from Eq. (1), see Fig.[3l Every crossing point gives 
the electronic temperature T e in the normal metal at a 
particular bias. Let us note that this comparison can- 
not be made on the differential conductance curve since 
there would be a systematic error in the determination 
of T e (V) due to its bias dependence. 

Fig. 2] displays the effective electronic temperature in 
the normal metal for Sample 2. Taking into account a 
leakage resistance of down to 5 MQ (which is the mini- 
mum possible as based on the I-V characteristic measured 
at 200 mK), would not change the extracted electronic 
temperature T e (V) by more than 1 mK at the gap edge. 
As bias is increased, the electronic temperature decreases 
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FIG. 3: (Color online) Sample 2 experimental current-voltage 
characteristic at a 275 mK cryostat temperature (full line, 
black) superposed on a series of calculated isotherm charac- 
teristics following Eq. (1) (dotted lines) from T = 292 mK 
(top) to T = 97 mK (bottom) . Every crossing point gives the 
electronic temperature T e in the normal metal at a particular 
bias. 



FIG. 4: (Color online) Sample 2 dependence of T e (dots) on 
the cooler bias for a cryostat temperature of 275, 470 and 570 
mK. The base temperature extracted from the data is respec- 
tively 292, 489 and 586 mK. The coloured full and dotted 
lines give the calculated electron and phonon temperatures 
respectively obtained from the electron-phonon model for E 
= 2 nW./im- 3 .K- 5 , KA = 66 pW.K~ 4 . 



until the bias hits the gap. This minimum temperature 
of about 100 mK for a base temperature of about 300 mK 
is typical for electronic coolers probed by conventional N- 
I-S thermometers [3[. At a bias above the gap, the elec- 
tronic temperature cannot be accurately determined the 
same way since the current- voltage characteristics calcu- 
lated at different temperatures overlap strongly. 

Fig. [5] inset shows the schematic of the thermal model 
we consider, which is similar to the one described in 
in the framework of hot electrons in metals. The normal 
metal electrons are at thermal equilibrium at the temper- 
ature T e and cooled by the two superconducting reser- 
voirs thermalized at the base temperature Tb ase . The 
heat current P coo i out of the normal metal into an indi- 
vidual N-I-S junction is given by: 

Pcooi = -5^- I" (E-eV/2)n s (E)[f N (E-eV/2)-f s (E)}dE 

6 tin J — oo 

(2) 

where fs is the energy distribution function in the super- 
conductor. The electron-phonon scattering is the domi- 
nant energy exchange driving energy into the electronic 
system with a power: 

P e _ ph (T e ,T ph ) = SC/(T e 5 - T* h ) (3) 

where E is a material dependent prefactor and U is the 
volume of normal island. Eq. [3] is derived under the hy- 
pothesis that the phonons have an energy spectrum sim- 
ilar to the bulk regime and that their energy distribution 
is thermal at a temperature T p h. Its validity was recently 
demonstrated experimentally in Cu submicron wires [12[ . 
According to [13], the actual distribution function will 



be out of equilibrium, with a lower energy fraction corre- 
sponding to the base temperature Tb ase and a higher en- 
ergy fraction at the electron temperature T e . We verified 
that this out-of-equilibrium distribution function can be 
approximated with an equilibrium one with an effective 
phonon temperature T p h in between Tb ase and T e . The 
heat balance equation for the N-metal electrons reads 
2Pcooi + Pe-ph — 0, where the factor 2 accounts for the 
presence of two symmetric N-I-S junctions. A power dis- 
sipation term due to a leakage resistance (of down to 5 
MVt which gives a power of only 2 fW) could be added to 
the above equation, but would not make any difference 
in the following discussion. 

The power transfer between the phonons of the normal 
metal and the substrate is given by the Kapitza coupling: 

P K (T ph , T base ) = KA(Tp h - T b \ se ), (4) 

where A is the interface area and K depends on the ma- 
terials in contact [14| . Here we consider that the phonons 
of the substrate always remain at the base temperature 
Tb ase . The heat balance equation for the normal metal 
phonons is given by P e - P h + Pk = 0. 

It is generally accepted that the phonons of a metal film 
are strongly thermalized to the substrate, corresponding 
to the limit of large K in the above model, leading to 
T p h = Tb ase . In this case, Eq. [3] implies that the quantity 
(T e /T hase ) b depends linearly on the quantity 2P coo i/T^ ase 
with a slope equal to 1/TXJ. Fig. [5] shows that although 
the low bias data at a given base temperature follows rea- 
sonably such a behaviour, the related values of the pa- 
rameter E depend significantly on the base temperature. 
Moreover, they are in the range 0.8-1.2 nW./im -3 .K -5 , 
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FIG. 5: (Color online) Sample 2 measured bias dependence 
of the quantity (T e /Tbase) 5 for a base temperature Tbase of 
292, 489 and 586 mK from bottom to top as a function of 
P /Tbase- Linear fits are displayed which gives values for the 
E parameter equal to 1.22, 1.02 and 0.78 respectively. Inset 
shows the schematic of the thermal model. 



i.e. well below the expected value of 2 nW./im _3 .K -5 
IE 0, EH • Thus we cannot fit our data within the ther- 
malized phonon hypothesis. Let us point out that, in 
contrast with previous studies of hot electron effects, we 
consider here a regime with a relatively high bath tem- 
perature and electrons cooled below that temperature. 

In order to describe our experimental results, we pro- 
pose to take into account a phonon cooling. We have 
self-consistently numerically solved the two heat balance 
equations as a function of sample bias in order to calcu- 
late T e (V) and T p h(V) for every base temperature. In 
the following, we have taken the electron-phonon cou- 
pling constant E to be equal to the well-accepted value 
E = 2 nW./im" 3 .K -5 and used the quantity KA as the 
fit parameter. We also used the measured values of the 
superconducting gap A and of the normal state resistance 
R n . We obtained a very good fit of the bias dependence 
of the electron temperature at every base temperature 
for KA = 66 pW.K~ 4 (Sample 2), see Fig.U At the gap 
edge, it is less accurate, presumably because of mech- 
anisms that are not included in our simple model, in- 
cluding the back-flow of hot quasi-particles accumulated 
close to the junctions. Similar results were obtained at 
every base temperature on all of our cooler devices with 
different tunnel resistances and junction areas. 

Fig. [H also shows, for three different base temperatures, 
the calculated effective phonon temperature as a function 
of the sample bias. At a base temperature of 489 mK, 
the phonons reach at the optimum bias (with a cooling 
power of 1.5 pW) a temperature about 50 mK below. At 
lower temperature, the ratio between the phonon cooling 
and the electron cooling decreases because the electron- 



phonon decoupling (oc T -4 ) dominates the Kapitza resis- 
tance (ex T -3 ). Above the gap, the electron temperature 
and the phonon temperature increase, since in this regime 
the tunneling brings heat to the electrons. As a consis- 
tency check, we calculated the differential conductance 
(Eq. (1)) of every cooler as a function of the bias, us- 
ing the bias-dependent electron temperature determined 
above. Fig. [2] shows Sample 1 calculated differential con- 
ductance curve compared to the experimental data. The 
two curves show a very good agreement (note the loga- 
rithmic scale). 

The parameter KA accounts for the total Kapitza cou- 
pling of the Cu island to the Si substrate, that includes 
the contributions of both a direct Cu/Si contact and a 
contact through the oxidized Al film with a presumably 
weaker efficiency. Dividing the fit-derived KA parame- 
ter by the full area of the Cu island provides an effective 
mean value (K). We found that (K) varies from sample 
to sample in the range 40 - 60 W.m -2 .K -4 . It is about 
a factor 3 below the values calculated for bulk Si/Cu in- 
terface [l4| and only a factor 2 below the experimental 
value in Ref. [y. 

In conclusion, the precise investigation of S-I-N-I-S 
coolers characteristics led us to study cold electron ef- 
fects with a relatively high temperature bath, in the 
framework of models previously mainly used to describe 
hot electrons effects with a very low temperature bath. 
We demonstrated the relevance of the phonon cooling, 
which can be characterized by the dimensionless param- 
eter YidT /K here of the order of 1. The description of the 
relevant thermal couplings between the phonon and elec- 
tron populations at the nanometer scale calls for further 
experimental and theoretical investigation. 

The authors are grateful to Nanofab-CNRS for sam- 
ple fabrication, E. Favre-Nicollin for contributing to 
early stages of this work and to P. Brosse for techni- 
cal support. We acknowledge useful discussions with J. 
Pekola, P. Gandit, F. Giazotto and I. J. Maasilta. This 
work was funded by EU STREP project 'SFINx' and by 
NanoSciERA project 'Nanofridge'. 



[1] M. Nahum, T. M. Eiles and J. M. Martinis, Appl. Phys. 

Lett. 65, 3123 (1994). 
[2] F. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin 

and J. P. Pekola, Rev. Mod. Phys. 78, 217 (2006). 
[3] M. M. Leivo, J. P. Pekola and D. V. Averin, Appl. Phys. 

Lett. 68, 1996 (1996). 
[4] A. M. Clark, N. A. Miller, A. Williams, S. T. Ruggiero, 

G. C. Hilton, L. R. Vale, J. A. Beall, K. D. Irwin and J. 

N. Ullom, Appl. Phys. Lett. 86, 173508 (2005). 
[5] Z. Jiang, H. Lim, J. Eom and V. Chandrasekhar, Appl. 

Phys. Lett. 83, 2190 (2003). 
[6] M. L. Roukes, M. R. Freeman, R. S. Germain, R. C. 

Richardson and M. B. Ketchen, Phys. Rev. Lett. 55, 422 

(1985). 



5 



[7] F. C. Wellstood, C. Urbina and J. Clarke, Phys. Rev. B 

49, 5942 (1994). 
[8] J. P. Pekola, D. V. Anghel, T. I. Suppula, J. K. 

Suoknuuti, A. J. Manninen and M. Manninen, Appl. 

Phys. Lett. 76, 2782 (2000). 
[9] D. R. Heslinga and T. M. Klapwijk, Phys. Rev. B 47, 

5157 (1993). 

[10] J. P. Pekola, T. T. Heikkila, A. M. Savin, J. T. Flyktman, 
F. Giazotto and F. W. J. Hekking, Phys. Rev. Lett. 92, 



056804 (2004). 

[11] H. Pothier, S. Gueron, D. Esteve and M. H. Devoret, 

Phys. Rev. Lett. 73, 2488 (1994). 
[12] J. T. Karvonen, L. J. Taskinen and I. J. Maasilta, J. of 

Low Temp. Phys. 146, 213 (2007). 
[13] N. Perrin and H. Budd, Phys. Rev. Lett. 28, 1701 (1972). 
[14] E. T. Swartz and R. O. Pohl, Rev. Mod. Phys. 61, 605 

(1989). 



